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Abstract 

Based upon our previous works ( Eur.Phys.J.D 6, 319(1999); Chin. Phys. Lett. 
18, 236(2001)), we develop a classical approach to calculate the high-order 
harmonic generation of the laser driven atoms and molecules. The Coulomb 
singularities in the system have been removed by a regularization procedure. 
Action-angle variables have been used to generate the initial microcanonical 
distribution which satisfies the inversion symmetry of the system. The nu- 
merical simulation show, within a proper laser intensity, a harmonic plateau 
with only odd harmonics appears. At higher intensities, the spectra become 
noisier because of the existence of chaos. With further increase in laser in- 
tensity, ionization takes place, and the high-order harmonics disappear. Thus 
chaos introduces noise in the spectra, and ionization suppresses the harmonic 
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generation, with the onset of the ionization follows the onset of chaos. 
PACS numbers: 42.65.Ky, 32.80.Wr, 32.80.Rm. 
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I. INTRODUCTION 



The development of the high-power femtosecond laser has stimulated the investigation of 
the multi-photon processes of atoms and molecules interacting with intense laser fields [1-9] . 
Recently, there are many theoretical and experimental references about these multiphoton 
process. Within a proper intensity region, lots of odd harmonics of laser are generated by 
atomic and molecular gases [1-5,14]. The harmonic structure distributes as a plateau which is 
cut off at a special high-order harmonic. When the intensity increases, the ionization channel 
is opened and the high-order harmonics disappear. The above threshold ionization (ATI) 
occurs when the laser intensity is sufficient power (above lO^^W / cm'^) . Various structures 
(plateau, angular distribution, etc.) in ATI spectra have been detailed in Ref.[6-8]. The 
above threshold dissociation (ATD) and dissociation-ionization of molecular systems are 
also reported [9-10]. Studying the laser-matter interaction deeply, not only can obtain new 
knowledge of the interacting mechanism, but also can provide widely application in the 
generation of high-order coherent harmonics, X-rays laser and 7-rays laser. 

The classical dynamics of most laser-driven systems is generally chaotic, due to the ex- 
istence of nonlinearity. Chaos usually manifests itself as some control parameters (initial 
energy, laser intensity, laser frequency, etc.) are varied. The microscopic systems, in partic- 
ular those involving atoms and molecules, are governed by a Hamiltonian. To study these 
systems are of great importance in the context of quantum-classical correspondence. At the 
investigated high power laser intensity {10^^~ 10^^ W/cm'^), the electric field of the laser is 
equal in strength to the Coulomb field of the nuclei [8] , and the laser field can not be looked 
as a perturbation to the field-free system, for the states themselves are no longer indepen- 
dent of the laser. Generally, the microscopic systems (atoms or molecules) are intrinsically 
quantum mechanical systems, thus they must be described by quantum mechanics. How- 
ever, due to the presence of intense laser, the exact calculation even numerical simulation 
based upon quantum mechanics tend to be very difficult to perform. Fortunately, classical 
approach to these similar problems is useful for providing physical insight into dynamics 



3 



processes [10-13]. Classical chaos associating with the microwave ionization of atomic hy- 
drogen reveals that the detailed mechanism of atomic ionization in terms of transport in 
phase space [11-12]. Classical prediction for scaled frequencies has been verified by quan- 
tum calculations [12-13], and, in turn, has confirmed the dynamical significance of classical 
chaos. The harmonic generation (HG) of laser driven hydrogen atoms is simulated with 
classical Monte-Carlo method [14], the numerical results are qualitatively consistent with 
the quantum mechanic results and the experimental observation. For the molecular systems, 
the classical calculation is a good first step, since length and energy scales are often large 
enough for classical mechanics to be at least approximately valid. Using softened potential 
model, the classical dynamics of the one-dimensional hydrogen molecular ion H2 interacting 
with an intense laser pulse are detailed [15-16]. 

There are two important and realistic examples of laser driven systems, where nonlinear 
dynamics has played a major role, which have been traditionally studied. One is the laser 
driven hydrogen atom that is the fundamental system in atomic physics, the other is the laser 
driven hydrogen molecular ion H2, which is the fundamental diatomic molecular system in 
molecular physics. Bellow, we shall show our classical calculation of the high-order harmonic 
generation of these two fundamental systems. The outline of this paper is as follows. The 
regularized model and the corresponding initial microcanonical distribution are presented in 
the next section. In section III, we show the numerical results in details. A briefiy summary 
and discussion is given out in the last section. 

II. MODEL AND INITIAL MICROCANONICAL DISTRIBUTION 

A. Regularized model of laser driven hydrogen atom 

In this article, we consider a classical hydrogen atom, with an infinite mass nucleus fixed 
at the origin of coordinates, interacting with a high intense laser field which is linearly po- 
larized along the z — axis, and with the electric field component s{t). Thus the Hamiltonian 
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in atomic units in Cartesian coordinates is 

H = Ho + Hi^ 

Ho^\p'-\.Hi^-ze{t). (1) 

where, r = \/x^ -\-y^ + z'^ and — pI + + pI- Apparently, the above Hamiltonian has a 
Coulomb singularity corresponding to electron-nucleus collision. To remedy the singularity, 
we introduce the parabolic coordinates {u, v, (p) 

X — uv cos (l),y — uv sin (f),z — {v"^ — ii^)/2. (2) 

and a new Active time scale r 

dt/dT^G{u,v)^v'^ + u'^. (3) 

In order to implement regularization, following the notation of Szebenhely [17], regarding the 
motion time t and the negative total energy —E as generalized coordinate and generalized 
momentum respectively, then the Hamiltonian function in extended space becomes into 

H* ^ H - E = ^ H*(x,y,z,t;p,,py,p,,-E) (4) 

For the sake of obtaining the regularized Hamiltonian, we introduce the third category 

generating function F^{pxiPyiPz-i —E; u, v, (f), t), then obtain 

dFs dFs dF3 dF^ 

x{u, V, 0) = —^,y{u, V, 0) = —g^^ ^(«' ^) = ^(^' ^) ^ ~ ■ 

So the third category generating function F2{px,Py,Pz, ~E;u,v,(l),t) is in the form of 

F3 = -xpx - ypy - zpz + Et, 

— —uv cos (ppx — uv sin (ppy — -(f^ — u^)Pz + Ft. (6) 

The new momenta ipu,Pv,P<j), —E) are related with old momenta ipx,Py,Pz, —E) by 

Pu = = V COS (pPx + V Sm (j)Py - upz, 

ou 



dv 
dFs 



Pv — E — ~ ^'^^ ^Px + 4'Py ~ '^Pz-i 



P4> = — = -uv sin (t)px + uv cos 
E^?§^E. (7) 



Then the regularized Hamiltonian can be expressed as 

K = ^H* = G(u, v)(H -E) = 0, 
dr 

-1\pI + pI + iu-' + v-')pI] - 2 - E{u' + v') -^-{v'- u')e{t). (8) 

If we define 

Ku = IIpI + u-y,] -1-Eu' + lu'e{t), 

K.^l\pl + v-y,]-l-Ev'^hMt). (9) 

In the field-free case, e{t) — 0, Ky,{— —K^) is the component of the Laplace-Runge-Lenz 
vector along z — axis, i.e., — —Ky — R^. The equations of motion can be derived from 
the above regularized Hamiltonian, as the following: 
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dK 
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dr d(f) 
2Eu+pIu-^ -2u^6{t), (11) 
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(12) 



and 



dE dK 1 4 4.de{t) dt dK 

-;— — — — v) — ; — , — = ———= G{u,v). (13) 
dr dt 2^ ^ dt ' dr dE ^ ' ^ ^ ^ 

Apparently, is a constant which corresponding to the component of angular momentum 
along z — axis. 



B. Regularized model of laser driven hydrogen molecular ion H2 

Within Born-Oppenheimer approximation, we can assume that two protons are fixed at 

the positions A and B, with a distance R away from each other, in the locations (0, 0, ~R/2) 
and {0,0, R/2) of the Cartesian coordinates system. A single electron at {x,y,z) is subject 
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to the Coulomb attraction of both protons and the interaction of the laser. Let ta be its 
distance from A, and tb be its distance from B. Within the dipole approximation, the 
Hamiltonian in atomic units is 

H = Hq + 

Ho = ^p' - l/rA - I/VB, Hi = -ze{t). (14) 

where, e{t) is the electric field of the laser pulse. Hi is the interacting Hamiltonian. Ap- 
parently, the above Hamiltonian has singularities at points — and = 0, which 
corresponds to electron-proton collision. To overcome this barrier, regularization has to be 
performed. Define the new coordinates {u, v, (f)) as 

R R R 

X — —— sinii sinh V cos 4>,y = —— sinw sinh vsin0, ^ = cos w cosh v. (15) 

and introduce the new Active time scale r satisfying 

dt/dr = G{u, v) = r^re = i?^(cosh^ v - cos^ u) /4. (16) 

Similar to the previous subsection, regarding the motion time t and the negative total energy 
—E as generalized coordinate and generalized momentum respectively, then we can write 
the Hamiltonian function in extended space as 

H* = H-E = = H*{x,y,z,t;p,,py,p„ -E) (17) 

Introducing the third category generating function F^{pxiPyiPz-i —E; u, v, (p, t), which satisfies 

R R R 

F3 = — — sin u sinh v cos (pp^ sinu sinh v sin (ppy cos u cosh vp^ + Et. (18) 

2 2 2 

Thus momenta ipu,Pv,P4>, ~E) are related with old momenta ipx,Py,Pz, ~E) by 

Pu — — -Q— = cos u sinh v cos (f)px + ^ cos m sinh v sin — ^ sin u sinh vpz , 
Pv — — -Q— = —smu cosh V cos 0Pa; + —smu cosh v sm (ppy — —smu cosh vp^ , 

R . • . • ^ ^ ^ • • . ^ 

Prj) — — -^T" — —-^z-smu smh v sm (pp^ + — smu smh v cos (ppy, 

E=^=E. (19) 



And the regularized Hamiltonian is 

K = ^H* = G{u, v)(H -E) = 0, 
dr 

= ■^\Pu + pI + (csc^ u + csc/i^vjp^] — i^coshv — ^i?^(cosh^ v — cos^ u). (20) 

where, Hi = —zs{t) = {—R/2)cosucosh.ve{t). Obviously, is a constant, which corre- 
sponds to the component of the angular momentum which along the z — axis. If is equal 
to zero, the electron is constrained in a plane which can be chosen as y — 0, corresponding 
to a two-dimensional motion. In the field- free case, e{t) = 0, this two-dimensional model 
corresponds to the classical hydrogen molecular ion in ground-state. Equations of the mo- 
tion for the ground-state hydrogen molecular ion interacting with laser field can be easily 
obtained from the regularized Hamiltonian, i.e., 

du _dK _ 
dr dpu 

dpu dK EE? . ,^ . , dG ^dz. 
dv dK 



dr dpy 

^ = = i^sinh^; + — — sinh 2^; +£ i 0^— + G— , 22 

dr ov 4 ov av 



and 



dE dK deit) dt dK 
d7 = ^ = d^^-dE^ 

In the field-free case, defining 

(pI + pI ^) + co^^ ^' 

1 EE? 
Ky = - (pI + p'^ CSC h'^v) — R cosh v — cosh^ v. (24) 

thus Ku{= —Ky) is a constant of the field-free motion, they are related to 7 and fl by 

X„ = -Ky = -7i?V4 = ER^A + n/2. (25) 



Constants Q and 7 are first introduced by Erickson [18] and Strand [19] respectively, which 
have the following forms 

J = -E -2n/{mR^), 

fl — La ■ Lb + emF?{cos9A — cosOb)- (26) 

Here, La and Lb are the angular momentum vector of the motion around nucleus A and 
B respectively, 9a and 9b are the angle from the vector rX and fs to positive z — axis 
respectively, m and e are the mass and the charge of the electron respectively. 



C. Action-angle variables and initial distributions 

The action-angle variables for separated systems are defined as 

1 / , . (9 /■ , r dpi dH , 
h-^f Pdq,, Oi=gj-J Pidq. - J gj^gj:dq,. (27) 

The motion of the classical field-free hydrogen atom is periodic, then it need only a pair of 

conjugated action-angle variables 



Ho = Eo = (28) 

The corresponding angle is given by the Kepler's equation 

9 = u — esinw. (29) 

This means 9 is the mean anomaly of the free orbits. The eccentric anomaly u can be 
obtained from r = a(l — ecosu). Here, r is the distance from the electron to the origin, a is 
the instantaneous semimajor axis which satisfies a — — l/(2£'o), e is the eccentricity of the 



free orbits satisfying e = ^J2EqL'^ + 1, L is the total angular momentum. 

The motion of the ground-state hydrogen molecular ion H2 is quasi-periodic, it need two 
pairs of action-angle variables. With eUiptic integrals, action 7„ can be expressed as follows. 
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for Ku > 0, 



TT 



u 



for = 0, 



(30) 



V-2goR2 



)], /or X„ < 0. 



TT 



where, the first category elhptic integral Fi{ip,k) and the second category elliptic integral 
F2{(fi, k) are in the form of 



Generally, a single trajectory lacks the inversion symmetry of the real physical sys- 
tems. So the spectrum obtained from a single trajectory exhibits unphysical even har- 
monics. A natural way to remedy the unphysical even harmonics is to consider an en- 
semble of trajectories, evolving from an initial microcanonical distribution with inver- 
sion symmetry. For a chaotic system, the initial distribution can be generated with 
Monte-Carlo method. However, for an integrable system, the distribution generated by 
Monte-Carlo method does not possess of ergodicity. We find that the points on the 
same equienergy surface generated by action- angle variables with regular steps possess 
of good ergodicity. To reconstruct the inversion symmetry, there must exist pairs of 
{xo,yo,ZQ;pci:o,Pyo,PzQ) and {-xq, -yo, -zq] -p^^, -py^, -p^o) in the initial distribution. For 
the regularized model for hydrogen atom, it corresponds to pairs of {uo,vo,4>o;PuotPvojP(I>o) 
and {uo, fo, 0o + '^'■iPuo:Pvo:P(i>o)- And for the regularized hydrogen molecular ion i?2^, it 
corresponds to pairs of {uo,Vo,(j)o;Puo,Pvo,P4>o) and {uq + ^:,Vo,(f)o;puo,Pvo,P<t>o)■ 



The numerical computational procedure is based upon the classical trajectory Monte 
Carlo (CTMC) method [20-21]. CTMC simulation procedure involves three stages, (i) choice 
of initial conditions, (ii) numerical integration of equation of motion, and (iii) categorization 













III. NUMERICAL SIMULATION 
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of each trajectory as excitation, charge transfer or ionization. In the process of numerical 
integrating, numerical accuracy and computing time are two primary aspects that must be 
considered. We use the fourth-order Runge-Kutta method with variable steps to perform 
the numerical calculation. Note also that computer can not deal with singularity, which 
corresponds to electron-nucleus collision, to overcome this difficulty, we have implemented 
regularization. With complete regularization, numerical simulation can be established with 
required precision before, at, and, after collision successfully. 

To obtain the harmonic spectra, our procedure also calculate the averaged dipole moment 
of the excited trajectories with pairs of inversion symmetric initial conditions in the same 
distribution. Having determined the actual trajectories of the electron, one can easily obtain 
the component of the averaged dipole moment, which along the laser polarization direction. 
Then the harmonic spectra of the driven dipole is straightforwardly obtained from it's power 
spectra 



D(uj) — lim - 







(31) 



where, (//(t)) is the component of the averaged dipole moment along the laser polarization 
direction. In our calculation, the electric fields of the ultrashort laser pulses are chosen as 



EMsm'^{uLt/4:0) sm{u;Lt), for 0<t< AQ-k/ujl, 

(32) 

0, otherwise. 



e{t) 



The maximum electric field strength Em is related to the laser intensity / = ^e^^j \xqE\iI2^ 
and lol\& the angular frequency, and period is equal to 2t: /ujl. 

The time-dependent dipole moment //^(t) of a single trajectory along the laser polariza- 
tion direction sensitively depends on the laser intensity. Fig.l shows the time evolution of 
the dipole of the hydrogen molecular ion with initial energy Eq — —1.1034 hatree, inter- 
nuclear distance it! = 2.00 bohr, laser wavelength A = 600 nm and different laser intensity. 
And some dipoles of the hydrogen atom and their power spectra arc presented in Fig.2, 
with initial energy Eq — —0.5 hatree, laser wavelength A = 532 nm and different laser 
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intensity. With the increasing of the laser intensity, the oscillations of the dipole moments 
are modulated gradually, and it follows a regular pattern with both high and low frequency 
components. Such patterns are independent of initial conditions. When the laser intensity is 
large enough, ionization takes place, which strongly modifies the subsequent time evolution 
of the dipole moment. 

The power spectra obtained from a single trajectory of the laser driven hydrogen molecu- 
lar ion, which with initial energy Eq — —1.1034 hatree, internuclear distance it! = 2.00 bohr, 
laser wavelength A = 600 nm and different laser intensity, are presented in Fig. 3. The first 
one corresponds to the field-free case. For the hydrogen atom, it possesses of only regularly 
decreasing peaks locating at multiples (harmonics) of Kepler's frequency, that is, cUn — ncuo, 
for n = 1, 2, 3 • ••, which manifests that the free motion is periodic. While for the hydrogen 
molecular ion, the free motion is quasi-periodic, it appears regularly decreasing peaks locat- 
ing at two different characteristic frequencies and their combinations, i.e., at ncuoi + TnujQ2-i 
for n = 0, 1, 2, 3 • ••, m = 0, 1, 2, 3 • •• and n -|- m > 0. For the laser driven systems, the peaks 
of harmonic spectra depend strongly on the laser intensity. In addition to the origin peaks, 
a dominant line at the laser frequency ujl appears, which is the Rayleigh component in the 
light scattered by the atoms and molecular ions. High-order harmonics, which consist of 
both odd and even components, do appear in the spectra, and their orders and strengths 
increase with the increasing of the laser intensity. For a low laser intensity, i.e., in the 
perturbative regime, the characteristic peaks dominate the power spectrum. For a proper 
laser intensity, the plateau structure containing both odd and even harmonics appears. For 
a larger laser intensity, because of the presence of chaos, the spectra become noisier even 
the motion is bound. For a strong enough laser intensity, the motion becomes unbound and 
the corresponding spectrum is dominated by a very noise background, which is generated 
by the ionized electrons, the only line surviving being the one at the laser frequency, which 
corresponds to the hght scattered by the asymptotically free electron Thomson scattering. 

To eliminate the unphysical even harmonics, averaging an ensemble of trajectories evolv- 
ing from an inversion symmetric distribution is necessary. The time evolution of the averaged 
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dipole, which evolve from a microcanonical ensemble of 5000 trajectories, are presented in 
Fig. 4, the first one corresponds to the hydrogen atom with initial energy Eq = —0.5 hatree, 
laser wavelength A = 532 nm and laser intensity / = 5.0x10^'' W/ cw? , the second one corre- 
sponds to the hydrogen molecular ion with initial energy Eq — —1.1034 hatree, internuclear 
distance R — 2.00 hohr, laser wavelength A — 600 nm and laser intensity / — 1.0 x 10^^ 
W/cm^. Comparing with the evolution of the dipole of a single trajectory, one can easily 
find that the oscillation of the averaged dipole is smooth and it globally follows the laser 
oscillation. 

The spectra obtained from an ensemble of trajectories are showed in Fig. 5 and Fig.6. 
Fig. 5 is the harmonic spectra of the ground-state hydrogen atoms interacting with the laser 
pulses with A = 532 nm and different laser intensity. Fig.6 is harmonic spectra of the 
hydrogen molecular ion with initial energy Eq = —1.1034 hatree, internuclear distance 
it! = 2.00 bohr, laser wavelength A = 600 nm and different laser intensity. As a consequence 
of averaging process, the unphysical even harmonics are remedied really. Within a proper 
laser intensity range, the plateau structure that only possesses of odd harmonics appears. As 
pointed out previously, at a higher intensity, the spectra become noisier even the ionization 
does not happen because of the effects of chaos. This indicates that the chaos cause the 
noise of the harmonic spectra. When the laser intensity is high enough, the ionization takes 
place, thus the noise background conceals the high-order harmonics. This means that the 
onset of ionization follows the onset of chaos and the ionization suppresses the harmonic 
generation. 

IV. SUMMARY AND DISCUSSION 

In summary, within the Born-Oppenheimer approximation and using the classical trajec- 
tory method,we have calculated the high-order harmonic generation spectra of the hydrogen 
atom and the hydrogen molecular ion interacting with ultrashort intense laser pulses. The 
other multi-photon phenomena, such as multi-photon ionization and above threshold disso- 
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ciation, can also be simulated. In our subsequent calculations, the dynamics of the electron 
is investigated by numerical integrating the equations of motion using regularized coordi- 
nates. To eliminate the unphysical even harmonics of a single trajectory, averaging over an 
ensemble of trajectories evolving from an initial microcanonical distribution with inversion 
symmetry is necessary. Such distribution is constructed using action-angle variables. A 
plateau structure in the spectra with only odd harmonics is observed within a proper laser 
intensity range of about 10^^ W/cm^. Prom our numerical results, we observe that the high 
order harmonics are cut off at a special order harmonic. At higher laser intensities, chaos 
introduces noise into the spectra even though the motion is still bound. Finally as the 
intensity is further increased, ionization takes place, and the harmonics disappear. 

These results are qualitatively consistent with recent quantum calculations [16] and ex- 
perimental observations [1-5], but the cutoff order A^^ of the plateau structure is not precisely 
consistent with formula — {Ip + 3.17Up)/hu!^[27], here, Ip is the ionization potential and 
Up — e^Efj/AmeCU^ denotes the quiver energy or the ponderomotive energy of an electron. 
As an example, when the laser intensity / = lO^'' W/cm^ and wavelength A = 600 nm, 
the ionization potential Ip of the ground-state hydrogen molecular ion is 1.1034 hatree 
(29.77ey = 14.50 fiwj, the quiver energy Up = 3.36eV (= 1.63 hwj, then A^^ = 19, and 
when / = 7.5 x 10^^ W/cm"^, Nm = 18; however, from our simulation, the harmonic plateau 
are both cut off at 17, and for very large intensity (above 10^^ W/cm'^) the high-order har- 
monics are concealed by the noises. To obtain quantitative results, we have to integrate the 
time-dependent Schrddinger equation, this can be realized with the split-operator method 
[22]. For the hydrogen molecular ion, within BOA, a proper internuclear distance R will 
enhance the high-order harmonic generation [16], and it will be interesting to go beyond the 
Born-Oppenheimer approximation to investigate what further interesting insights can be 
obtained when the nuclear motion is taken into account [23-24]. 

In our model, we only consider the non-relativistic case with dipole approximation. When 
the laser is sufficiently intense, photoelectrons of relativistic energies can be produced, ne- 
cessitating a fully relativistic treatment [25]. The dipole approximation is no longer valid, 
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and the magnetic field is not only present, but acquires an importance similar to that of the 
electric field. And even before the appearance of the relativistic photoelectrons, the effects 
of the magnetic field may be very important too [26]. Due to the wiggly motion and the 
acceleration of the electron near the outermost turning points induced by the magnetic field, 
the cut-off order of the harmonic plateau maybe higher. 
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FIGURES 

FIG. 1. Temporal variation of the dipole moment of a single electronic trajectory of the hy- 
drogen molecular ion for different laser parameters with initial energy Eq = —1.1034 hatree and 
internuclear distance R = 2.00 bohr. 

FIG. 2. The dipole and their power spectra of a single electronic trajectory of the hydrogen 
atom for different laser parameters with initial energy Eq = —0.5 hatree. 

FIG. 3. Power spectra of the dipole of a single electronic trajectory of the hydrogen molecular 
ion for different laser parameters with initial energy Eq = —1.1034 hatree and internuclear distance 
R = 2.00 bohr. 

FIG. 4. Temporal variation of the averaged dipole moment for different laser parameters. This 
first one corresponds to the hydrogen atom with initial energy Eq = —0.5 hatree, the other 
one corresponds to the hydrogen molecular ion with initial energy Eq = —1.1034 hatree and 
internuclear distance R = 2.00 bohr. 

FIG. 5. Harmonic spectra of the ground-state hydrogen atoms interacting with different laser 
pulses, the last two are magnifications of the first one. 

FIG. 6. Harmonic spectra of the ground-state hydrogen molecular ion interacting with different 
laser pulses, with initial energy Eq = —1.1034 hatree and nuclear internuclear distance R = 2.00 
bohr. 
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